Ising model in small-world networks 
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The Ising model in small- world networks generated from two- and three-dimensional regular lattices 
has been studied. Monte Carlo simulations were carried out to characterize the ferromagnetic tran- 
sition appearing in these systems. In the thermodynamic limit, the phase transition has a mean-field 
character for any finite value of the rewiring probability p, which measures the disorder strength of 
a given network. For small values of p, both the transition temperature and critical energy change 
with p as a power law. In the limit p — > 0, the heat capacity at the transition temperature diverges 
logarithmically in two-dimensional (2D) networks and as a power law in 3D. 
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I. INTRODUCTION 

Complex networks describe many systems in nature 
and society, and have been modeled traditionally by ran- 
dom graphs [jO. In the last years, new models of com- 
plex networks have been introduced, motivated by em- 
pirical data in different fields 0. In particular, small- 
world networks have been studied, as they are suitable to 
describe properties of physical systems with underlying 
networks ranging from regular lattices to random graphs, 
by changing a single parameter g]. Watts and Strogatz 
j| proposed a model for this kind of networks, which 
is based on a locally highly connected regular lattice, in 
which a fraction p of the links between nearest-neighbor 
sites are randomly replaced by new random links, thus 
creating long-range "shortcuts" . The networks so gen- 
erated are suitable to study different kinds of physical 
systems, such as neural networks Q and man-made com- 
munication and transportation systems ^jH^- 

These small-world networks interpolate between the 
two limiting cases of regular lattices {p — 0) and random 
graphs (p = 1). In the small- world regime, a local neigh- 
borhood is preserved (as for regular lattices), and at the 
same time some global properties of random graphs are 
maintained. The small-world effect is usually measured 
by the scaling behavior of the characteristic path length 
£, defined as the average of the shortest distance between 
two sites. For a random network one has a logarithmic 
increase of £ with the network size TV (i.e., the number 
of sites) , while for a e?-dimensional regular lattice one ex- 
pects an algebraic increase: £ ~ N x / d . In contrast, for a 
small- world network, £ follows the scaling law |)| |ll| 



£(N,p) ~ (N*) 1/d F(N/N*) 



(1) 



p 1 is a 



where the scaling function F(u) has the limits F(u 
u 1/d for u < 1 and F(u) ~ lnu for u > 1; TV 
crossover size that separates the large- and small-world 
regimes |],|l^,[l3|] . This indicates that the small- world 



behavior appears for any finite value of p (0 < p < 1) 
as soon as the network is large enough, and in particu- 
lar the global characteristics of the network are changed 
dramatically in the presence of only a small fraction of 
shortcuts. 

The importance of this shorter global length scale has 
been studied for several statistical physical problems on 
small-world networks. Among these problems, one finds 
in the literature the signal propagation speed the 
spread of infections 0,^5) , as well as site and bond per- 
colation Jlf]-[l7j. The localization-delocalization transi- 
tion of electron states has been also studied on quantum 
small- world networks [H . 

Up to now, most of the published work on small 
worlds has concentrated on networks obtained from one- 
dimensional lattices. Barrat and Weight PJ and Gitter- 
man |l9f have studied the crossover from ID to mean- 
field behavior for the ferromagnetic Ising model, which 
presents a phase transition of mean-field type for any 
value of the rewiring probability p > 0, provided that 
the system size is large enough. Close to p — the tran- 
sition temperature T c goes to zero as l/\ logp|. A mean- 
field-type behavior has been also found for the XY model 
in small-world networks generated from one-dimensional 
chains pOjj . 

More recently, Svenson and Johnston |MJ have studied 
the damage spreading for Ising models on small-world 
networks obtained by rewiring two-dimensional (2D) and 
three-dimensional (3D) regular lattices. They found that 
these networks are more suitable than regular lattices to 
study social systems with the Ising model. 

Here we investigate the ferromagnetic transition for 
the Ising model in small-world networks generated by 
rewiring 2D and 3D lattices. Contrary to the networks 
generated from ID lattices, now one has phase transitions 
at finite temperatures T c > for p = 0. This means that 
one expects a change from an Ising-type transition at 
p = to a mean-field-type one in the small-world regime. 
We employ Monte Carlo (MC) simulations to obtain av- 
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FIG. 1. Fourth-order Binder's cumulant Un as a function 
of temperature for small-world networks generated from 3D 
cubic lattices with rewiring probability p = 0.02. Symbols 
represent different system sizes N — L 3 : squares, L = 15; 
triangles, L = 20; circles, L — 25; diamonds, L = 30. 

erage magnitudes for finite-size systems. The resulting 
quantities are extrapolated to the thermodynamic (in- 
finite size) limit, where the small- world behavior is ex- 
pected to dominate the thermodynamic properties for 
any finite value of the rewiring probability p > 0. 

In Sec. II we describe the computational method. In 
Sec. Ill we present results of the MC simulations along 
with a discussion. The paper closes with some concluding 
remarks in Sec. IV. 



II. COMPUTATIONAL METHOD 



We consider the Hamiltonian 

H = ^ JijSiS 



3 > 



(2) 



where Si = ±1 (i = 1, ...,N), and the coupling matrix 
Jij is given by 



J i ■ 



_ J J(> 0), if i and j are connected, 



0. 



otherwise. 



(3) 



Monte Carlo simulations have been carried out on net- 
works of different sizes, generated from the 2D square 
and 3D cubic lattices. Small-world networks were gener- 
ated by randomly replacing a fraction p of the links of 
the regular lattices with new random connections. This 
procedure keeps constant the total number of links in 
the rewired networks. Thus, the average coordination 
number z in the 2D and 3D cases amounts to 4 and 6, 
respectively. In the rewiring process we avoided isolated 
sites (with zero links). With this procedure we obtained 
networks in which more than 99.9% of the sites were con- 
nected in a single component (Note that a random graph 
has usually many components of various sizes). The re- 
maining sites (when they appeared) were excluded from 
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FIG. 2. Transition temperature T c , normalized by the 
average coordination number 2, as a function of the rewiring 
probability p for small- world networks generated from 2D and 
3D lattices. Dashed lines are guides to the eye. 



the final networks employed for the MC simulations. The 
size of the networks used in our calculations was larger 
than the crossover size N* Jl]|[l(| , so that we were in the 
small- world regime. 

The largest networks employed here included 200 x 
200 sites for the 2D system and 40 x 40 x 40 sites for the 3D 
network. Periodic boundary conditions were assumed. 
Sampling of the configuration space has been carried out 
by the Metropolis local update algorithm (22|. Several 
thermodynamic quantities and moments of the order pa- 
rameter have been calculated for different simulation-cell 
sizes. Finite-size scaling was then employed to obtain the 
magnitudes corresponding to the thermodynamic limit 
(extrapolation to infinite size). In the remainder of the 
paper, the presented values for the different quantities 
will correspond to the extrapolated ones, unless explicit 
mention is made indicating a particular network size. 

The ferromagnetic transition temperature has been de- 
termined by using Binder's fourth-order cumulant p2] 
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where the magnetization M of a given spin configuration 
is given by M = Si/N. As an example, we present 

in Fig. 1 the cumulant Un as a function of temperature 
for different system sizes in the 3D case for p = 0.02. 
The transition temperature is obtained from the unique 
crossing point for the different sizes N. 

The heat capacity per site, c v , was obtained from the 
energy fluctuations at a given temperature, by using the 
expression 



(AE) 2 



Nk B T 2 ' 



(5) 



where {AE) 2 = (E 2 ) - (E) 2 
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FIG. 3. Dependence upon the rewiring probability p of 
the shift in transition temperature AT C with respect to the 
regular lattices, for 2D and 3D networks. Lines are guides to 
the eye. 



III. RESULTS AND DISCUSSION 

In Fig. 2, we present the transition temperature T c 
as a function of the rewiring probability p for 2D and 
3D networks. Similar to the ID case [E0[, T c changes 
fast close to p = 0, and the derivative dT c /dp becomes 
smaller as p increases. However, the sharp change of T c 
for small p shows a behavior different than that found in 
the ID case, where T c ~ | logp| _1 . In the limit p = 1 one 
finds an increase in T c as z rises from 4 (in 2D) to 6 (in 
3D), as expected for random lattices, for which one has 
ksTc/zJ — ► 1 as z — ► oo. 

To analyze the change in critical temperature with p, 
we call AT C — T c — Tj?, being Tj? the transition temper- 
ature for the corresponding 2D or 3D regular lattice. In 
Fig. 3 we show the dependence of AT C upon p for the 2D 
and 3D cases in a log-log plot. In both cases we find that 
AT C can be fitted by a power law AT C ~ p s for p < 0.01. 
The exponent s is 0.52 ± 0.03 for 2D and 0.96 ± 0.04 
for 3D. Our result for 2D networks is compatible with a 
y/p dependence for AT C near p = 0, which means that 
the derivative dT c /dp diverges as ~ 1/^/p for p — > 0. In 
the 3D case, our results seem to indicate a dependence 
AT C ~ p for small p. However, in this case the lowest 
p values studied here may still be too high to attain the 
small-p regime (see below). 

Associated with the increase in T c as the rewiring prob- 
ability p rises, one expects an increase in the critical en- 
ergy E(T C ). We call e c = E(T C )/N the critical energy 
per site, and Ae c = e c — e° c its change with respect to the 
regular lattice (p = 0). This difference Ae c is shown in 
Fig. 4 as a function of p for 2D (squares) and 3D (cir- 
cles) networks, in a log-log plot. For p < 0.01, Ae c can 
be well fitted in both cases by a power law of the form 
Ae c ~ p u . For the exponent u, we find u = 0.43 ± 0.03 
and 0.56 ± 0.04 in 2D and 3D, respectively. 

A characterization of the ferromagnetic phase transi- 



FIG. 4. Dependence on the probability p of the shift in 
critical energy per site Ae c with respect to the regular lattices, 
for 2D and 3D networks. Lines are guides to the eye. 



tion in these networks requires the determination of the 
universality class to which it corresponds. In the limit 
p = (regular lattices), one has transitions of the 2D 
and 3D Ising type. In order to determine the type of the 
phase transition at p > 0, we have studied the critical 
exponent (3, which gives the temperature dependence of 
the order parameter close to the transition temperature: 
(M) ~ (T c - T)P for T < T c . For the different values of 
the rewiring probability p studied here, we have calcu- 
lated numerically the logarithmic derivative 



M(t) 



rflog(M) 
dlogt 



(6) 



for t = T c — T > 0, which is related to the exponent j3 
through the limit f3 — lim t ^o 

In Fig. 5 we present results for the derivative fi as a 
function of temperature for several values of p and for 
a 2D network of size 200 x 200. For reference, we also 
present results of MC simulations for p — (Ising model 
on a regular 2D lattice) for the same system size, which 
converge to (3 — 0.125, the critical exponent for the 2D 
Ising model. In all cases p > 0, the extrapolation T — > T c 
gives an exponent (3 close to 0.5, the value corresponding 
to a mcan-field-type transition. However, for decreas- 
ing p, the observation of the mean-field character of the 
transition requires to go to temperatures closer to T c (or 
equivalently to larger system sizes; see 0]). Thus, for 
p = 0.001 we are still far from the value 0.5 at the tem- 
peratures at which the employed system size allows us to 
give a precise value for the order parameter (M). But 
even in this case, the departure from the value expected 
for a 2D Ising-type transition is clear close to T c . 

A complementary way to confirm the mean-field char- 
acter of the phase transition consists in determining the 
exponent is, which controls the behavior of the correla- 
tion length near the critical temperature: 



(7) 
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FIG. 5. Logarithmic derivative [i, versus the temperature 
difference t = T c — T for small-world networks generated by 
rewiring a 2D lattice of size 200 x 200. Different symbols 
represent results obtained for several values of the rewiring 
probability p. From top to bottom: p — 1 , 0.1, 0.01, 0.001, 
and 0. 



Close to T c , this critical exponent is related to the tem- 
perature dependence of the fourth-order cumulant de- 
fined in Eq. (|) as MM: 



U N (T) ^U* + U 1 



1 fl " 



(8) 



with U* and U\ independent of T and the system size N. 
From this expression we have: 



At/ 



N 



AT 



oc 



(9) 



which allows us to calculate the exponent v from the cu- 
mulant Un derived from the MC simulations. The values 
of v obtained by this method agreed in all cases (within 
error bars) with the critical exponent corresponding to 
mean-field transitions: v = 0.5. 

We conclude that the mean-field character of the phase 
transition (which is the one found in random networks) 
should appear in the thermodynamic limit for any p > 
0. This is in line with the observation mentioned above 
that the characteristics of the random graphs (e.g., mean 
distance between sites I ~ log N) show up in small- world 
networks as soon as their size N is large enough. 

To understand the dependence on p of the transition 
temperature T c for small p, we will consider the two 
length scales present in this problem. On one side we 
have the correlation length £, which near T c follows the 
temperature dependence given in Eq. (Q). On the other 
side, we have a length scale characteristic of the small- 
world network, given by the typical distance between 
ends of shortcuts: £ = {p>z)~ x l d |l6). When the corre- 
lation length is smaller than £, the system behaves ba- 
sically as a regular lattice. When £ grows beyond £ (as 
happens when we approach the transition temperature 
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FIG. 6. Heat capacity per site c v versus temperature 
for small-world networks generated from a 2D lattice of size 
200 x 200. The plotted curves correspond to different values 
of the rewiring probability p, as indicated by the labels. 

from above, T — > T+), the "long-range" interactions in- 
troduced by the shortcuts come into play and give rise 
to the mean-field behavior. Thus, the transition between 
the regular-lattice behavior and the mean-field one oc- 
curs for £ ~ C- For the ID Ising model, the correla- 
tion length at low temperatures (kgT <C J) is given by 
£ ~ exp(2J/kBT) j24|. Taking into account that in this 
case £ ~ l/p, the condition ^ k ( suggests a critical 
temperature for the ID small-world T c ~ |logp| _1 , in 
agreement with more detailed calculations for this sys- 
tem For the 2D Ising model the correlation length 
near T c scales as £ - \T - T^- 1 HQ, whereas now 
C ~ l/p 1 / 2 . Then, using the same argument, one ex- 
pects for small-world networks generated from 2D lat- 



tices: T r - T > 



P 



1/2 



in good agreement with the re- 



sults of our Monte Carlo simulations (see Fig. 3). Thus, 
it seems that in general, for systems with a correlation 
length diverging as \T — T®\~ v , the order-disorder tran- 
sition temperature for small p depends on the rewiring 
probability as 



T r - T j 



p 



l/vd 



(10) 



Note that here v is the critical exponent of the considered 
model in the regular lattice, not the one corresponding to 
the phase transition in the small-world network (which is 
the mean-field one). Hence, we see that the condition £ « 
£ describes the transition between the large- world regime 
(which corresponds to regular lattices) and the small- 
world one 13 1, as well as the ferromagnetic transition, 
indicating that the small-world transition and the order- 
disorder one are intimately related. 

The above argument, however, seems to fail for small- 
world networks built up from 3D lattices. For the 3D 
Ising model one has a critical exponent v « 0.63 pi], 
which should give for small p a dependence: AT C ~ p ^ 6 , 
with an exponent clearly smaller than that derived from 
our MC simulations, where we found for p < 0.01: 



P 



0.96 



It seems that this discrepancy appears be- 



cause the smallest p value employed in our simulations 
(p = 10 -3 ) is still too large to observe the small-p be- 
havior. Although for such values of p we are clearly in 
the small-world regime, the corresponding value of C in 
3D (C ~ 5), is too small to allow the correlation length 
to reach the critical dependence of the Ising model in 
the regular lattice: £ ~ |T — T®\~ v . In other words, to 
observe the p-dependence of T c given in Eq. (|l0|), one 
needs £ values larger (i.e., smaller p values or larger net- 
works) than those employed here. A similar conclusion 
was proposed for the XY model on ID small-world net- 
works [^0| , for which MC simulations gave a dependence 
T c a logp+6 (with a and b numerical constants) instead 
of T c ~ p, expected from the above argument. 

For mean-field-type transitions, the heat capacity per 
site c v shows a finite jump (not a divergence) at T c . Given 
that at p — one has Ising-type phase transitions with 
divergent heat capacity at T c , such a jump in c v will 
diverge in the limit p — > 0. The temperature dependence 
of c v is displayed in Fig. 6 for several values of p and 
for networks built up from a 200 x 200 2D lattice. As 
expected, one observes an increase in the maximum value 
of c v as p is reduced. 

The maximum value of c v for each p, extrapolated to 
the thermodynamic limit, is shown in Fig. 7 for 2D and 
3D small-world networks. In Fig. 7(a) we present our 
results for 2D, which for small p display a logarithmic 
dependence c™ = a — bhip, with the numerical constants 
a = 1.74 ± 0.06 and b = 0.20 ± 0.01. This logarithmic 
dependence of the heat capacity can be explained by ar- 
guments similar to those employed above to analyze the 
behavior of the transition temperature as p — > 0. For 
the Ising model in 2D regular lattices, the heat capacity 

near T c ° diverges as: c v (T) log \T - T°| Taking 

into account the relation between T c and p given above 
for small p, one finds for the 2D small-world networks 
c™ ~ — logp, in agreement with the results of our Monte 
Carlo simulations. 

In the 3D Ising model, c v has at T? a singularity of 
the form |T - T?\- a , with a m 0.12 §1. Assuming that 



T r - T? 



P 



l/ud 



, as in Eq. (|10|), one expects for c™ close 



to p = a power-law of the form: c™ ~ p w , with an 
exponent w = —a/vd = —0.063. In Fig. 7(b) we present 
the values of c™ derived from our MC simulations as a 
function of p in a log-log plot. These results are consis- 
tent with a power-law dependence c™ ~ p w at small p. In 
fact, for p < 0.01 our results can be fitted with an expo- 
nent w = — 0.057 ± 0.005. However, as in the case of the 
transition temperature discussed above, smaller values of 
p are necessary to determine unambigously this depen- 
dence from numerical simulations, and in particular to 
find the exponent w. 
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FIG. 7. Maximum value of the heat capacity per site c™ as 
a function of the rewiring probability p. Symbols correspond 
to the extrapolation of finite-size results to N — > oo. (a) 2D 
networks, in a semi-logarithmic plot; (b) 3D networks, in a 
log-log plot. Dashed lines are guides to the eye. 
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IV. CONCLUDING REMARKS 

We have studied the ferromagnetic transition that ap- 
pears for the Ising model in small-world networks gener- 
ated from 2D and 3D regular lattices. In these networks, 
the presence of a small disorder (p > 0) causes a change 
in the universality class of the order-disorder transition, 
from Ising for p = to mean- held type for p > 0. 

Our results indicate that the order-disorder transition 
occurs at a temperature T c where the spin correlation 
length £ is on the order of the length £ (typical distance 
between ends of shortcuts), characteristic of these net- 
works. In particular, close to p — this gives a depen- 
dence T c ~ | logp^ 1 in ID networks, and T c -T c ° - p x l vd 
for networks generated from regular lattices of higher di- 
mensions. This is the dependence found from our Monte 
Carlo simulations for 2D networks. In the 3D case we 
find a power law for T c — Xj?, but the determination of 
the actual exponent from MC simulations for p — > re- 
quires p values smaller (i.e., networks larger) than those 
employed here. 

From the results presented in this paper, it is clear that 
there is a close relation between the small-world tran- 
sition and the order-disorder transition in this kind of 
networks. 

The author benefited from useful discussions with M. 
A. R. de Cara and M. Saboya. Thanks are due to 
E. Chacon for critically reading the manuscript. This 
work was supported by CICYT (Spain) under Contract 
No. PB96-0874, and by DGESIC through Project No. 
1FD97-1358. 
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